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We construct equilibrium configurations of uniformly rotating neutron stars for selected relativis¬ 
tic mean-field nuclear matter equations of state (EOS). We compute in particular the gravitational 
mass (M), equatorial (i?eq) and polar (Apoi) radii, eccentricity, angular momentum (J), moment of 
inertia (1) and quadrupole moment (M 2 ) of neutron stars stable against mass-shedding and secular 
axisymmetric instability. By constructing the constant frequency sequence / = 716 Hz of the fastest 
observed pulsar, PSR J1748-2446ad, and constraining it to be within the stability region, we obtain 
a lower mass bound for the pulsar, Mmin = [1.2-1.4] Mq, for the EOS employed. Moreover we give 
a fitting formula relating the baryonic mass (M;,) and gravitational mass of non-rotating neutron 
stars, Ms/Mq = M/Mq -I- (13/2OO)(M/M0)^ [or M/Mq = Ms/Mq — (1/2O)(M6/M0)^), which is 
independent on the EOS. We also obtain a fitting formula, although not EOS independent, relating 
the gravitational mass and the angular momentum of neutron stars along the secular axisymmetric 
instability line for each EOS. We compute the maximum value of the dimensionless angular mo¬ 
mentum, a/M = cJ/{GM^) (or “Kerr parameter”), (a/M)max ~ 0.7, found to be also independent 
on the EOS. We compare and contrast then the quadrupole moment of rotating neutron stars with 
the one predicted by the Kerr exterior solution for the same values of mass and angular momen¬ 
tum. Finally we show that, although the mass quadrupole moment of realistic neutron stars never 
reaches the Kerr value, the latter is closely approached from above at the maximum mass value, as 
physically expected from the no-hair theorem. In particular the stiffer the EOS is, the closer the 
Kerr solution is approached. 
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I. INTRODUCTION 

Understanding the physics of neutron stars and being 
able to describe the spacetime in the interior and around 
such compact objects is one of the most important objec¬ 
tives for modern astrophysics, as from one side it could 
be an excellent test for the general relativity and on the 
other because the strong forces acting and the high level 
of densities reached cannot be tested in any place in the 
earth. 

Although our actual incapability to describe these 
concepts exactly, some observational limits has already 
been determined, via simplifying assumptions (such as 
a spherical configuration in the X-ray binaries observa¬ 
tion). Other general considerations on the nature of neu¬ 
tron stars and pulsars are often extracted in the liter¬ 
ature from the use of fiducial structure parameters: a 
canonical neutron star (NS) of mass M = 1.4 Mq, radius 
i? = 10 km, and moment of inertia / = 10“^® g cm^ (see, 
e.g., [l|, and references therein). Based on these param¬ 
eters together with the so inferred surface magnetic field 
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from the classic point-like magneto-dipole rotating model 
[ii , NSs have been traditionally classified according to 
the thought nature of the energy source powering their 
observed emission. Neutron stars are also thought to 
possibly participate in the most powerful explosions in 
the Universe, the gamma-ray bursts (GRBs), e^^. via 
NS mergers in the case of short GRBs (see, e.g, Q, and 
references therein) and hypercritical accretion processes 
leading to gravitational collapse to a black hole (BH) in 
the case of long GRBs associated to supernovae (see, e.g., 
[ii, and references therein). 

There are still, however, many open issues regarding 
the above global picture of NSs both from the physics 
and astrophysics points of view (see, e.g., 0,i). On the 
other hand, our theoretical and observational knowledge 
on NSs has been largely increased in the intervening years 
from the first general relativistic description of a NS by 
Oppenheimer and Volkoff in 1939 Q. Eventually a more 
complex equation of state (EOS), interior structure, and 
consequently exterior gravitational field of non-rotating, 
slowly rotating, and fast rotating rotating stars have been 
acquired, and massive NSs of ~ 2 Mq, drastically con¬ 
straining the nuclear EOS stiffness, have been observed 
MM- Thus, general conclusions based on fiducial pa- 

observables, might be pre- 
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mature and a more exhausting exploration of the conse¬ 
quences of adopting different NS parameters appears to 
be necessary. 

For instance, as recently shown in [l^ for the high- 
magnetic field pulsar class, thought to be intermediate 
objects linking pulsars and magnetars, the understand¬ 
ing of their magnetic field values and their energy source 
needs the precise knowledge of the entire set of possi¬ 
ble structure parameters of a NS and a self-consistent 
general relativistic description of the pulsar observables. 
Indeed, it was there shown how the magnetic field of a 
pulsar is overestimated, and the rotational energy under¬ 
estimated, by the classic magnetic rotating dipole model 
and the use of hducial NS parameters. 

The description of the rotational and thermal evolu¬ 
tion, as well as the emitted radiation of isolated and ac¬ 
creting NS requires the knowledge of the structure prop¬ 
erties and the corresponding exterior metric. For in¬ 
stance, we have recently compared and contrasted the 
cooling evolution of neutral NSs satisfying local charge 
neutrality and global charge neutrality [l^. We have 
shown there that, owing to their different crust struc¬ 
ture (mass and thickness) for a same value of the total 
mass, their thermal relaxation time (the time to form an 
isothermal core) can be very different and therefore the 
signatures of the structure of the NS might be accessed 
via early cooling observations. 

There are recent numerical relativity’s computations 
of the structure of uniformly rotating NSs, which have 
mainly focused on the existence of universal relations 
between e.g. the quadrupole moment, moment of in- 
ertia,, and Love number of NSs (see, e.g., Bin , and 
[la considered in a slow rotation regime. Other 

works (e.g. Bill) tried to recover these relation also 
in a full rotation regime, with a numerical method based 
on the work firstly implemented in [^ . 

In this work, through a full rotation approach (treated 
by numerical relativity’s methods), we will focus on ad¬ 
ditional structure properties of uniformly rotating NSs 
relevant for astrophysical applications such as mass, po¬ 
lar and equatorial radii, eccentricity, angular momentum, 
angular velocity, moment of inertia, and quadrupole mo¬ 
ment, for a selected sample of EOS (describing nuclear 
matter with relations of different stiffness) within rela¬ 
tivistic mean-field nuclear theory, not analyzed in the set 
of EOS of previous works. 

This article is organized as follows. In section [IT] 
we briefly review the axisymmetric system of Einstein’s 
equations to be integrated for a given EOS which we 
describe in section Ell The stability conditions (mass- 
shedding and secular instability) are outlined in section 
m and the mass-radius of rotating NSs is shown in sec¬ 
tion El The eccentricity and the moment of inertia are 
shown in section Ell while the quadrupole moment is dis¬ 
cussed in section Em We finally summarize and discuss 
our results in section ElIIl 


II. STRUCTURE EQUATIONS 

We consider the equilibrium equations for a self- 
gravitating, rapidly rotating NS, within a fully general 
relativistic framework. We start with the stationary ax¬ 
isymmetric spacetime metric (see, e.g., [llj): 

ds^ = + e^'^{dcj) — Ludt)^ + e^^{dr^ + r^dO^), ( 1 ) 

where v, ^/;, w and A depend only on variables r and 0. 
It is useful to introduce the variable = rsm{d)Be~'^, 
being again B = B{r,9). The above form of the metric 
is obtained under two assumptions: 1) the Killing vec¬ 
tor fields are, one timelike and one relative to the 
axial symmetry, (/>“; 2) the spacetime is asymptotically 
flat. Then, one can introduce quasi-isotropic coordinates, 
which in the non-rotating limit they tend to isotropic 
ones. 

Turning to the physical matter content in the NS inte¬ 
rior, if one neglects sources of non-isotropic stresses, vis¬ 
cosity, and heat transport, then the energy-momentum 
tensor becomes the one of a perfect-fluid 

= (e-HP)u“n^-f P5“^, (2) 

where £ and P denote the energy density and pressure of 
the fluid, and u°‘ is the fluid 4-velocity. In terms of the 
two Killing vectors, 

V1 — V 

where v is the fluid 3-velocity with respect to the local 
zero angular momentum observer (ZAMO), 

v={n-uj)e^-' (4) 

being Q = u'^/u^ = d(j)/dt the angular velocity in the 
coordinate frame, equivalent to the one measured by an 
observer at rest at infinity. 

Thus, with the metric given by equation (HD and the 
energy-momentum tensor given by equation @ , one can 
write the field equations as (analogoulsly to Ref. [1^ set¬ 
ting C = A -I- i^): 


V • {BVv) = sin^ OB^e-'^^Vuj ■ Vw 


47ri?e^^ 


(e -I- P)(l -I- v"^) 
l-n2 


2P 


,( 5 ) 


V • (r^ sin^ 6»p3e“'‘'"Va;) = -167rr sin6»P' 

X (6) 

1 — 

V • (rsin(6»)VP) = 167rr sin6»Be^^"^'^P, (7) 
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- (l - + ^ (l - r‘^B~^B^rU},fj.U}^r - ^ (l - M^) r‘^B~^B^^j, r^iuj^rf - {^J?) {u},f,f | 


( 8 ) 


r 


where, in the equation for C./i, we introduced /r = cos(d). 

The projection of the conservation of the energy- 
momentum tensor, normal to the 4-velocity, 
((5g -I- u'^Ub) VaT°'^ = 0, leads to the hydrostation¬ 
ary equilibrium equation: 


1 




= 0 , 


Pi + {e + P) -f -- 2 ( -vv^i +v'^ 

^ ' 11 — LU / 

(9) 

where A = 1, 2, 3 and, as usual, A^i = dAjdx^. 

For a barotropic equation of state (EOS), P = P(e), 
and in the case of uniform rotation which we adopt in 
this work, the above hydrostationary equilibrium equa¬ 
tion has a first integral that can be written as 


dP 
e + P 


ln('U Vat) — i^lpole: 


( 10 ) 


where the constant of motion has been obtained, for in¬ 
stance, at the pole of the star (see, e.g., Ref. iH). 


III. EQUATION OF STATE 

To obtain a solution to the field equations, the matter 
EOS must be supplied. In general, a NS is composed of 
two regions, namely the core and the crust. The core, 
with densities overcoming the nuclear saturation value, 
Pnuc ~ 3 X 10^^ g cm“^, is composed by a degenerate gas 
of baryons (e.g. neutrons, protons, hyperons) and leptons 
(e.g. electrons and muons). The crust, in its outer region 
{p < Pdrip ~ 4.3 X 10^^ g cm“^), is composed of ions and 


electrons, while in its inner region (pd 


rip 


< P < Pn 


c), 


there is an additional component of free neutrons dripped 
out from nuclei. For the crust, we adopt the Baym- 
Pethick-Sutherland (BPS) EOS [2^. For the core, we 
here adopt modern models based on relativistic mean- 
field (RMF) theory. Indeed, RMF models have become 
the most used ones in NS literature, being its success 
mainly owing to important properties such as Lorentz 
covariance, intrinsic inclusion of spin, a simple mecha¬ 
nism of saturation for nuclear matter, and being con¬ 
sistently relativistic, they do not violate causality (see, 
e.g.. Ref. [15 ). We adopt, as now becoming traditional, 
an extension of the original formulation of Boguta and 
Bodmer in which nucleons interact via massive me¬ 
son mediators of different nature providing the attractive 
long range (scalar a) and repulsive short range (vector 
(jj) of the nuclear force, isospin and surface effects (vec¬ 
tor p). Meson-meson interactions can be also present; for 
instance, in the version of Boguta and Bodmer [2^ there 
is the presence of a self-interacting scalar field potential 
in form of a quartic polynom with adjustable coefficients. 
We consider here the possibility of including, in addition 
to such potential, vector-vector interactions of the w me¬ 
son. For a very recent and comprehensive analysis of the 
performance of several RMF models in the description of 
observed properties of ordinary nuclei, we refer the reader 
to Ref. [^, and for a brief historical and chronological 
reconstruction of the developments of the RMF models. 
Ref. 


27|. 


Thus, we shall constrain ourselves to models in which 
the energy density and pressure are given by (in units 
with ?i = c = 1) [ 2 ^: 


J 


1 2 2 I 9<y2 3 ffcr3 4 1 2 2 , 2 2\2 ^ 2 2 , , 9p . 

^ - -jJ^pPo + + yPo^s + 2^ 

i—n,p^e,^ 




79 1 2 2 ffo'2 3 9tT3 4,1 2 2 I 9(b3 , 2 2\2 ,1 22, 

P = -b —(ga;Wo) +^m^po-b 


E 

i—n,p,e,ii 


p, 


(lla) 

(llb) 


where 171 ^^ 01 ,p are the masses of the scalar and vector mesons, (?cr 2 , 3 i 9un 9 ui 3 are coupling constants, cr, Wq, po 
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denotes the scalar meson and the time-component of the 
w and p vector mesons, respectively. The components £i 
and Pi for each kind of particle considered are 


(27r)3 Jo 
2 

^ / 
/o 


k'^ + {m^p)'^cFk, (12a) 


£e,f, = j -k (me,^)2d^fc, (12b) 


P -i- 2 

n,p — - 


3 (27r)3 


0 + 


zdJk, (12c) 


Pe,ii = 


1 2 


3(27r)3_/g y^fc2 _|_ (rrie,^)' 


(12d) 


where with m* is intended the effective mass of 
baryons. 

The scalar, isospin, and baryon densities are given by, 
respectively. 


Us 


2 

(^ 



i—n,p 


m. 


\/fc2 + (to*) 


-.Pk, 


7 Z 3 — Hp Tin-, 

TIq — Tip Tlji. 


(13a) 

(13b) 

(13c) 


where, m = (fcf )^/(37r^), are the particle number densi¬ 
ties with kf the particle Fermi momenta. 

The equations of motion of the meson fields within the 
RMF approximation are: 


mla = gaUs - gcr2cr^ - gascr^, 
= 9LjnB - gui39ui{9uit^of, 

2 9p 
^pPO = Y«3- 


(14a) 

(14b) 

(14c) 


A barotropic EOS can be obtained iff additional clo¬ 
sure relations are supplied. A first condition to be im¬ 
posed is the request of the stability of matter against 
beta decay. The second closure equation that has been 
traditionally adopted is the condition of local charge neu¬ 
trality of the system. The latter condition has been re¬ 
cently shown to be not fully consistent with the equilib¬ 
rium equations in presence of multicomponent charged 
constituents such as protons and electrons (see [13 , and 
references therein). Instead, only global charge neutral¬ 
ity must be requested. The new system of equations, 
referred to as Einstein-Maxwell-Thomas-Fermi (EMTF) 
equations, introduce self-consistently the Coulomb inter¬ 
actions in addition to the strong, weak, and gravitational 
interactions within a full general relativity framework. It 
is worth to notice that in this case no perfect-like form 
of the total energy-momentum tensor is obtained since 
the presence of electromagnetic fields breaks the pressure 
isotropy. Static NSs fulfilling the EMTF equations were 
constructed in Ref. [13], and uniformly rotating configu¬ 
rations in the second-order Hartle approximation can be 
found in Ref. [13 ■ To construct rotating NSs beyond the 
slow rotation regime, we take advantage of existing public 


numerical codes (e.g. the RNS code, see section Ej) that 
solve the field equations without any limitation of the 
rotation rate of the star. However, an implementation 
of the equations and boundary conditions of the EMTF 
system within these codes is not yet available. Thus, as a 
first step, we adopt in here the condition of local charge 
neutrality, bearing in mind the necessity of a future im¬ 
plementation of the EMTF equations of equilibrium in 
fast rotation regime. 

With both the beta equilibrium and the local charge 
neutrality conditions, a numerical relation between the 
energy density and the pressure can be obtained. We 
here adopt the nuclear parameterizations (for the spe¬ 
cific values of the coupling constants, particle, and me¬ 
son field masses) NL3 [1|, TMl [13, and GMl 
(after having adapted these in units used by the code). 
In figure [I] we compare and contrast the three selected 
EOS used in this work in the nuclear and supranuclear 
regime, relevant for NS cores. Below in section El we 
shall show how this selection of EOS are not only physi¬ 
cally relevant, as it is shown in this section, but also from 
the astrophysical point of view. 



log(£/c^ [gem ^1) 


FIG. 1. Pressure-energy density relation for the three EOS 
(TMl,GMl, and NL3) used in the present article. 


IV. STABILITY OF EQUILIBRIUM MODELS 

To solve the system of field equations dSHS]) and the 
hydrostationary equilibrium equation m, one has to fix 
one (in the spherical static case) or two parameters (in 
the rotating case). The first quantity one has to fix is 
the central value of energy density, £c- For a rotating 
model, one can choose the second parameter amongst 
different possibilities: the axes ratio (rpoi/req) of coor¬ 
dinate radii, angular velocity (Q), dimensionless angular 
momentum (j), gravitational mass (M) or baryonic mass 
(M{,). Thus, it is always possible to construct a sequence 
of rotating models by fixing a value of the second param¬ 
eter and letting the central energy density to vary in a 
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given range which is constrained to stability limits which 
we now discuss. 

The hrst limit on the stability of uniformly rotating 
configurations which we take into account is given by the 
sequence of maximally rotating stars, also referred to as 
Keplerian or mass-shedding sequence. In all the stars be¬ 
longing to such a sequence, the gravitational force equals 
the centrifugal force at the star equator, in such a way 
that faster rotation rates would induce the expulsion of 
mass from the star. The RNS code calculate this se¬ 
quence by decreasing the axis ratio (which correspond to 
an increase in angular velocity) until the angular veloc¬ 
ity equals the one of a particle orbiting the star at its 
equator. 

Another limit to the physically relevant models is de¬ 
termined by the called secular axisymmetric instability. 
For static configurations, the maximum stable mass (the 
critical mass) coincides with the first maximum of a se¬ 
quence of configurations with increasing central density, 
namely the first point where dM/dSc = 0, being M the 
mass of a configuration with central density Eq- At this 
point, the frequency of the radial perturbations vanishes. 
For higher densities, imaginary frequencies are obtained 
which lead to non-oscillatory perturbations, hence an in¬ 
stability. Thus, for static configurations, a turning-point 
of the M-Ec relation locates the onset of unstable config¬ 
urations. This instability proceeds on secular timescales, 
i.e. not dynamical, so that it proceeds on long times that 
allow the star to accommodate itself to the energy loss 
that occurs when going from one equilibrium point to 
another during gravitational collapse (see, e.g., Ref. [s^ . 
and references therein). As shown by Friedman, Ipser 
and Sorkin in [s^, the turning-point method leading to 
points of secular instability can be also used in uniformly 
rotating stars as follows. In a constant angular momen¬ 
tum sequence, the turning point of a sequence of configu¬ 
rations with increasing central density separates secularly 
stable from secularly unstable configurations. Namely, 
secular axisymmetric instability sets in at 


dM (ec, J) 

OEc 


= 0 , 

J=constant 


(15) 


and therefore the curve connecting all the maxima 
(turning-points) limits the stability region. The intersec¬ 
tion of such a limiting curve with the Keplerian sequence 
gives the fastest possible configuration. It is important 
to mention that the numerical code adopted (described 
in the next section) builds sequences of constant dimen¬ 
sionless angular momentum, defined as 


■ _ 

^ ^ GM|’ 


(16) 


which is the quantity we refer to in the sequel. 

The angular momentum J is computed from the defi¬ 
nition 


J = [ 


(17) 


being E a spacelike 3-surface, = Vat/|V&tV^f| the 
unit normal vector field to the t=constant spacelike hy¬ 
persurfaces and dV = ^y\^g\d^x the proper 3-volume el¬ 
ement (with ^g the determinant of the 3-metric). With 
this, equation (HU becomes 

J = [ r^sin^{9)drd9d(t). (18) 

J l-v^ 


V. MASS-RADIUS RELATION, 
OBSERVATIONAL CONSTRAINTS, AND 
STABILITY REGION 

There are in the literature many different numerical 
schemes and consequently codes to compute relativis¬ 
tic, rotating figures of equilibrium. For the numerical 
integration of the equilibrium equations we use in this 
work the public code RN^ by Stergiuolas and Fried¬ 
man . This code is a numerical implementation based 
on the scheme by Cook, Shapiro and Teukolsky 
(firstly implemented for realistic NS EOS in [ 2 ^), which 
is a modified version of the method envisaged by Ko- 
matsu, Eriguchi and Hachisu [s^ . We refer the reader to 
Ref. for further details on the numerical schemes. 

The major intuitive effect of rotation is to deform the 
figure of equilibrium with respect to the spherical static 
counterpart. This can be seen from many points of view. 
Eor instance, we can compute sequences of constant an¬ 
gular velocity O. An important aspect should be taken 
into account nowever: the RNS code builds fast rotat¬ 
ing models starting from a spherical (static) guess and 
decreases the polar to equatorial radii ratio until the pa¬ 
rameter fixed (e.g. the angular velocity) is reached with a 
prescribed accuracy. Thus, the axes ratio is a parameter 
used intrinsically by the numerical method, while other 
parameters (see the beginning of previous section for a 
list) can be chosen, but are reached spanning decreas¬ 
ing values of axis ratio. In particular, as an example, 
the code does not converge for every value of fixed an¬ 
gular velocity in every range of central energy density, 
and the range of convergence gets reduced by decreas¬ 
ing the angular velocity. To be more precise, choosing 
fixed rotation frequencies below 300 Hz, the code fails to 
converge in the entire range in which equilibrium mod¬ 
els should exist (thus between the Keplerian and secular 
instability limits), even adopting a very dense numerical 
grid (300 angular times 600 radial points), different ac¬ 
curacy and tolerance values (down to lO”^^) or values of 
a relaxation factor from I to 0.8. Effectively the code 
does converge for this kind of rotation frequencies, but 
in very limited ranges of energy density. Thus, how can 
the slow rotation regime be recovered in all the stability 
region? As a technical advise, we would like to mention 
that one could compute first, sequences keeping various 


^ http://www.gravity.phys.uwm.edu/rns/ 
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values of axes ratio constant (in the vicinity of unity), 
and then select in this set of models the ones with small 
values of a particular angular velocity. We construct with 
this simple method the sequences with low rotation fre¬ 
quency (e.g. models from 50 Hz to 200 Hz). For other 
values of rotation frequency and for all other parameters 
constant sequences, we achieved optimal convergence us¬ 
ing a 300 angular times 600 radial points numerical grid, 
and accuracy and tolerance of 10“®, while the relaxation 
parameter was not necessary. 

Figure [2] shows the total mass-central energy density 
plane for the selected EOS TMl, GMl and NL3. We 
also show the stability limits discussed above in section 
llVI and show explicitly some constant angular momentum 
sequences. 

Figure |3] shows instead the total mass-central energy 
density plane but in this case we show explicitly some 
selected constant rotation frequency sequences ranging 
from 50 Hz all the way up to the rotation frequency of 
the fastest observed pulsar, PSR J1748 — 2446ad, with 
/ = H/(27r) « 716 Hz [13. 

In figure m we plot the same H—constant sequences to 
show the relation between M and the equatorial radius, 
77eq- 

With the knowledge of the mass-radius relations pre¬ 
dicted by the theory, we are now in the position of com¬ 
paring and contrasting them with existing observational 
constraints, in order to validate the selection of EOS of 
the present work, not only from the already presented 
physical aspects, but also from the astrophysical ones. 
Current observational constraints to the mass-radius re¬ 
lation of NSs are (see figure [5|): 

• Most massive NS observed. The mass value of 
the most massive NS observed is the one of PSR 
J1614-2230 with 2.01 ± 0.04Mq [s^- The rotation 
frequency of this pulsar is 46 Hz; thus the devia¬ 
tions from spherical symmetric are negligible. This 
implies that every mass-radius relation for non¬ 
rotating NSs must have a maximum stable mass 
larger than this value. 

• Fastest observed NS. The highest rotation fre¬ 
quency observed in a pulsar is the one of PSR 
J1748 - 2446ad with / = 716 Hz [H. The con¬ 
stant frequency sequence of this value for any mass- 
radius relation must have at least one stable con¬ 
figuration that supports such a rotation frequency; 
namely the constant frequency sequence for this 
pulsar must lie within the region of stability. This 
is actually a very weak constraint since most of NS 
models allow much higher rotation frequencies. In¬ 
terestingly, as we shall show below, the construc¬ 
tion of the constant frequency sequence for PSR 
J1748 — 2446ad allows us to infer a lower mass for 
this pulsar. 

• Constraints to the NS radius. Since the surface 
temperatures of not-so-young NSs {t > 10^-10^ y) 


are of the order of million degrees (see, e.g., [il), 
their thermal spectrum is expected to peak in the 
soft X-rays. Thus, the modeling of the NS X- 
ray emission appears to be, at the present, one 
of the most promising methods to obtain infor¬ 
mation on the NS radius. Systems which are 
currently used to this aim are isolated NSs, qui¬ 
escent low-mass X-ray binaries (qLMXBs), NS 
bursters, and rotation-powered millisecond pulsars 
(see. Ref. [13, and references therein). From the 
modeling of the observed spectrum, the radius of 
the NS as measured by an observer at infinity, 
i?oo = R/v'l - 2GM/(c2R), can be extracteco. 
The observation of a preferable radius at infinity 
clearly represents a constraint to the NS mass- 
radius relation since the above definition for R^o 
can be rewritten as 2GM/= R — R^/{R^). In 
Ref. im (see, also. Ref. [13), the X-ray emission 
from the NSs in the qLMXBs M87, NGC 6397, 
MI3, w Gen, and NGC 6304 was revisited, and 
in Ref. [13 the one of the NS X7 in the Globu¬ 
lar Cluster 47 Tucanae. From the extracted val¬ 
ues of Roo consistent with these observational data 
at 90% confidence level, we can conclude that the 
current X-ray data constraints very weakly the 
mass-radius relation, allowing radii in the interval 
Roo = [7.64,18.86], where the lower limit is ob¬ 
tained for NGC 6304 and the upper one for X7. It 
is important to mention that X-ray measurements 
suffer from a variety of uncertainties which are the 
main reason for the very large spread in possible 
NS radii. The spectra modeling depends on the 
atmosphere composition, magnetic fields, accurate 
knowledge of the distance to the source, hence the 
extinction, and to some extent on the NS exterior 
geometry which could be affected by the rotation 
of the NS in the case of some LMXBs which could 
harbor NSs rotating with frequencies of a few hun¬ 
dreds of Hz (see, e.g.. Ref. [13, for details). In these 
latter cases, a more reliable comparison between 
theory and the above data constraints, which as¬ 
sume spherical symetry, could be obtained by plot- 
ing mass-radius relation using, instead of the equa¬ 
torial radius, a mean or average spherical radius 
such as the authalic radius, {R) = (2i?eq + Rpoi)/3. 
However, for the purposes of this work, it is suffi¬ 
cient to make a comparison with the mass-radius 
relation produced by the non-rotating configura¬ 
tions. 

An additional constraint to the mass-radius relation 
might come from the request of causality to the EOS, 


^ Actually, accurate spectra modeling leads to preferable values for 
both mass and radius; however, for a simpler comparison between 
different results from different methods and for a simple test of 
the mass-radius relation it is sufficient to plot the constraints 
obtained from the values of i?oo consistent with the data M- 




7 


- Static sequence 

- Keplerian sequence 

- Secular Instability 



14.9 15.0 15.1 

log(ec/'^ [e 



J.^.» lU.U iu.i- 

log(ec/=^ [gcm^*]) 


log(ec/c^ [gem *]) 


FIG. 2. Gravitational mass is plotted against central energy density for j-constant seqnences obtained with the EOS TMl, 
GMl and NL3 (from top to bottom). In this plot and hereafter, the red, green and black curves represent respectively the static 
seqnence, the Keplerian sequence, and the limit for secular stability. Here other colors stand for various j-constant sequences. 


namely the condition that the speed of sound in the NS 
interior cannot exceed the speed of light. However, for 
the present set of EOS this condition is automatically 
satisfied by construction since the models are relativis¬ 
tic. One can therefore see from figure[5]that the spherical 
(static) models (solid curves) obtained by the EOS selec¬ 
tion of this work are well in agreement with the current 
constraints of the NS mass-radius relation determined by 
most updated observational data. 

An often useful physical quantity to be computed is the 
binding energy of the configurations, or the relation be¬ 
tween the baryonic mass and the gravitational mass. Eor 
non-rotating NSs, we found that for the three analyzed 
EOS, the following relations hold: 


Mb ^ 

M 

13 

1 

/ M 

Mq ' 

Mq 

200 


M 

. Ml 

1 / 

' Mb \ 


Mq 

20 1 

Mq) 


where Mh is the baryonic mass, hold, and thus appear 
to be a universal property. The maximum relative errors 
obtained for non-rotating sequences of GMl, TMl and 
NL3 are respectively 1.4%, 1.3% and 0.99%. Eor rotating 
configurations, M = M{Mt,J) or Mb = we 

find that for the our set of EOS there is indeed a common 
relation given by 


u - 1 . 7 ^ 

Mq Mq ^ 200 \Mq) 130-^ ) ’ 


( 20 ) 


accurate within an error of 2%, and which duly general¬ 
izes equation (HU). 

Turning back to the above plots, we can clearly see 
that, as expected, the higher the frequency of rotation 
is, the higher results the value of the mass at which start 
the departures from the non-rotating mass-radius rela¬ 
tion. We find that for rotation frequencies ^ 200 Hz (or 
rotation periods > 5 ms), the non-rotating star becomes 
an accurate representation of the object. This is in accor¬ 
dance with previous results, see e.g. figure 5 in Ref. [4l| . 
where it was shown that the moment of inertia of se¬ 
quences computed with different EOS, start to deviate 


considerably from the static and the slow-rotation Har- 
tle’s approximations for frequencies above ^ 0.2 kHz. As 
we show below, this is also the case for the moment of in¬ 
ertia in the same range of frequencies (thus, the moment 
of inertia of non-rotating configurations, can be safely ap¬ 
proximated with the one of spinning configurations, with 
frequencies below the aforementioned limit, and vicev- 
ersa). Eor higher frequencies, full rotation effects are 
needed for an accurate description. This is especially 
important for objects with masses lower that the maxi¬ 
mum value, where departures from a non-rotating or slow 
rotation approximation become more and more evident. 

Following this reasoning, it is important to see how 
a constant frequency sequence imposes structure con¬ 
straints to a pulsar. Particularly interesting becomes 
the case of the / = 716 Hz sequence (blue curve), 
which corresponds to the fastest observed pulsar, PSR 
J1748 — 2446ad. The constant frequency sequence inter¬ 
sects the stability region in two points: at the maximally 
rotating Keplerian sequence, defining a minimum mass 
for the pulsar, and at the secular axisymmetric insta¬ 
bility limit, in the upper part, defining the maximum 
possible mass for the given frequency. Clearly these min¬ 
imum and maximum mass values depend upon the EOS. 
For the EOS employed here, we can see that the mass of 
PSR J1748- 2446ad has to be > [1.41,1.35,1.20] Mq for 
NL3, TMl, and GMl, respectively. 

We now determine the maximum rotation frequency 
of NSs. The fastest configuration for a given EOS is the 
one that terminates the Keplerian sequence, namely the 
configuration at the intersection between the Keplerian 
and the secular axisymmetric instability sequences. We 
show in figure |6| the rotation frequency of the maximally 
rotating configurations, i.e. the frequencies of the NSs 
along the Keplerian sequence. 

Another important quantity for this discussion is the 
dimensionless angular momentum (“Kerr parameter”), 
a/M = cJ/{GM^), which we show in figure 0 as a 
function of the total mass for the maximally rotating 
configurations, namely the Keplerian sequence. It can 
be seen how the maximum value attained by the NS, 
(a/M)max ~ 0.7, holds for all the selected EOS. The max- 






























14.4 14.6 14.8 15.0 15.2 



log(£c/c^ [gem ^1) 

FIG. 3. Mass versus central energy density using the EOS 
TMl, GMl and NL3 (from top to bottom) set of parame¬ 
ters. Red and green curves represent the static and Keplerian 
sequences. Other colors correspond to constant frequency se¬ 
quences of value 716 Hz (fastest observed pulsar; blue), 500 Hz 
(purple), 300 Hz (orange), 200 Hz (red diamonds) and 50 Hz 
(black triangles). 


imum value is reached for the mass [0.96,1.05,3.37] Mq 
for the TMl, GMl, and NL3 EOS, respectively. The ex¬ 
istence of such a particular maximum, EOS-independent, 
value of a possibly implies the existence of universal lim¬ 
iting values of the NS compactness and the rotational 
to gravitational energy ratio. This is a conjecture which 
deserves further exploration. In the same plot, the same 
sequences obtained with other already known EOS (rep¬ 
resented by differently dashed curves), obtained assum¬ 
ing widely different kind of interactions and via differ¬ 
ent many-body-theories, are shown, and the reader can 



-Rcq [km] 




14 16 18 20 22 24 26 


[km] 

FIG. 4. Mass versus equatorial radius using EOS TMl, GMl 
and NL3 (from top to bottom) set of parameters for same 
sequences as in figure ([3]). The convention of the plot colors 
and symbols is the same as in figure (3) 


notice a general universal behavior of the dimensionless 
angular momentum, even if for these other EOS the ex¬ 
act maximal values of this dimensionless parameter are 
slightly different. On the other hand, such a general be¬ 
havior of a parameter in not surprising, in fact, as it was 
already shown in Ref. pjgj j , it can be chosen as a parame¬ 
ter to establish a universal I-Love-Q relation. Neverthe¬ 
less the important argument here is that, although differ¬ 
ent stiffness, our chosen set of EOS presents a common 
maximal dimensionless angular momentum a/M. 

In Table U we summarize a few relevant quantities of 
NSs; the maximum stable mass in the non-rotating case, 
the maximum mass in the case of uniform rotation, the 
maximum rotation frequency, and the maximum value of 
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FIG. 5. Mass versus equatorial radius relation using the TMl, 
GMl and NL3 EOS (respectively with colors red, blue and 
green), plotted together with up-to-date observational con¬ 
straints. Solid color curves represent static NS configurations; 
dotted color curves represent the sequence of models rotat¬ 
ing with spin frequency of the fastest observed pulsar (PSR 
J1748-2446ad), / = 716 Hz, and dashed color curves repre¬ 
sent sequences of models rotating at the Keplerian frequency. 
The gray-shaded region corresponds to the constraints given 
by the X-ray data while the horizontal lines are the lower and 
upper bounds to the mass of the most massive observed pul¬ 
sar, PSR J1614-2230, namely M = 2.01 TO.OdM©. Further 
details of these constraints can be found in the text. 



log(e/c^ [gem ^]) 

FIG. 6. Frequency of the maximally rotating configurations 
(Keplerian sequence) as a function of the total NS mass for 
the TMl, GMl and NL3 EOS. The curves end at the max¬ 
imum frequency configuration, which is located at the inter¬ 
section between the Keplerian and the secular axisymmetric 
instability sequences. 


the dimensionless angular momentum. 

Before closing this section, we would like to provide 
a formula, useful for astrophysical applications, for the 
masses of the NSs lying along the secular axisymmetric 
instability line. Using the dimensionless angular momen¬ 
tum j, defined in equation (nsi) and related to the Kerr 



FIG. 7. Dimensionless angular momentum (“Kerr parame¬ 
ter”), a/M = cJ/(GM^), as a function of the total NS mass 
along the Keplerian sequence for the EOS selected in this 
work (colored curves). For comparison we show the results 
for additional EOS, taken from a set supplied by the RNS 
code (EOS.INDEX file). We refer the reader to the RNS web 
page and references therein for further details of these EOS. 

TABLE I. Some properties of NSs for the selected EOS: crit¬ 
ical mass for non-rotating case, maximum mass in 

uniform rotation, maximum rotation frequency, /max, 

and maximum dimensionless angular momentum (“Kerr pa¬ 
rameter”), (o/M)max = [cJ/(GM^)]max. 


EOS 

[Mq] 

[Mq] 

/max [kHz] 

(a/M)max 

TMl 

2.20 

2.62 

1.34 

0.70 

GMl 

2.39 

2.84 

1.49 

0.71 

NL3 

2.81 

3.38 

1.40 

0.71 


parameter by j = {M/Mq)'^ a, we obtained 

M = Mi=°{l + kj^), (21) 

where the values of are given in Table HI A: = 

[0.017,0.011,0.0060] and I = [1.61,1.69,1.68] for the EOS 
TMl, GMl, NL3, respectively. The maximum relative 
errors obtained for values of mass along the secular ax¬ 
isymmetric instability line with respect to fits for each 
EOS are respectively [0.33%, 0.44%, 0.45%]. 


VI. ECCENTRICITY AND MOMENT OF 
INERTIA 

In order to see how a figure of equilibrium becomes 
deformed by rapid rotation, we compute the eccentricity 



which we plot as a function of the mass, M, for the same 
constant 0 sequences of the previous figures. 
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FIG. 8. Eccentricity versus gravitational mass using EOS with TMl, GMl and NL3 (from left to right) set of parameters for 
same sequences as in figure [S] 


It is also interesting to investigate the distribution of 
the energy density within the figure of equilibrium both 
the static and rotational case for the different EOS. In 
figure ini we show the contours of constant energy density 
of a model with central value Ec = 10^® g cm“^ both in 
the static case and in the rotational one with dimension¬ 
less angular momentum j = 4, for the sake of example 
for the GMI EOS. 

We turn now to compute the moment of inertia of the 
star, which is one of the most relevant properties in pulsar 
analysis. The moment of inertia can be estimated as [2l| 


where J is the star angular momentum which is given by 
equation (fT^ . 

In figure[Tn]we plot the moment of inertia as a function 
of the mass for some fl-constant sequences together with 
the Keplerian sequence, while in figure [TT] we show the 
relations between I and the compactness, 


where M* and i?* are the mass and the radius of the 
spherical configuration with the same central density as 
the rotating one, Ec- 

The above hgures conhrm that for rotation frequencies 
< 200 Hz, or rotation periods ^ 5 ms, the deformation of 
the star is very little and, indeed, the non-rotating or the 
slow rotation regimes can be safely adopted as accurate 
approximations of the rotating NS. 

VII. QUADRUPOLE MOMENT 

The quadrupole moment in the RNS code is given by 

= (WFi, (24) 

where req is the value of the coordinate radius at equator, 
p = 2v — ln(i?), s = r/(r -|- Teq) £ [0,1] is a compacted 
radial coordinate, p = cos(0), P 2 {p) is the Legendre poly¬ 
nomial of second order, and Sp = r'^Sp, being Sp a source 
function defined as 


Sp{r,p)= e2 


87re"^(£ + P) 


1 -I- 

1-U2 


,2„-2p 




+ Ge 


1\ 1 


w, 


+ ^(1 “ 


-7.r 






(25) 


with 7 = ln(i3). However, as shown in Ref. 0 , equa¬ 
tion (IM)) is not the actual quadrupole moment of the 
rotating source according to the Geroch-Hansen multi¬ 
pole moments (dsl-l^. Indeed, the quadrupole moment 


extracted via the Ryan’s expansion method [i^ is 00 


- M^l 


s'^ds' 
(I - s')® 


(26) 


(27) 

Jo 
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FIG. 9. Contours of constant energy density of a model with central value Ec = 10^® g cm ® both in the static case (left plot) 
and in the rotational one with dimensionless angular momentum j = 4 (right plot) for the GMl EOS. 





FIG. 10. Moment of inertia versus mass relation using EOS TMl, GMl, and NL3 (from top to bottom) set of parameters for 
same sequences as in figure [3] 





FIG. 11. Moment of Inertia versus compactness using EOS with TMl, GMl, and NL3 (from left to right) set of parameters 
for same sequences as in figure |3l 


where M 2 is given by equation (1^ and Tq is the Gegen- 
bauer polynomial of order 0 with normalization Tg' = 

^/2JttCo, with Co the traditional Oth-order Gegenbauer 
polynomial. 

Following Refs. [a [13, we computed numerically the 
correcting factor bo given by equation (HZD, and then ob¬ 
tained the corrected quadrupole moment through equa¬ 
tion (l26l) . In figure fT^ the modulus of is plotted 

in logarithmic scale against the gravitational mass for se¬ 
lected constant frequency sequences. Each sequence was 
stopped at the secular instability limit. We can see that 


the quadrupole moment is a decreasing function of the 
mass along a constant frequency sequence while it is an 
increasing function along the Keplerian sequence. 

We turn to compare and contrast the above mass 
quadrupole moment with the one of the Kerr solution. 




M' 


( 28 ) 


The reason for this is twofold. First, we point out the 
longly discussed question in astrophysics if the Kerr so¬ 
lution may describe the exterior gravitational field of a 
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- Keplerian sequence 

- /=300 Hz 

- /=500 Hz 

- /=716Hz 



FIG. 12. The modulus of corrected value for the mass quadrupole (in logarithmic scale) obtained via equation (I26II is plotted 
against gravitational mass for the same constant frequency sequences of figure [S] 


realistic astrophysical source besides a black hole, namely 
is there any matter content which could generate a Kerr 
exterior field? (see, e.g.. Refs, [i^, H^l, and references 
therein). Second, if the answer to the previous ques¬ 
tion is negative, then one can distinguish a NS from a 
black hole with the same mass and angular momentum 
from the knowledge of the quadrupole moment (see, e.g., 
Ref. [m, and references therein). 

In figure [13] we show the ratio between the NS 
quadrupole moment, given by equation dUl), and 

the Kerr solution quadrupole moment, for se¬ 

lected constant frequency sequences. We find that 
starts to approach as intuitively expected, for 

masses close to the maximum stable value. An inter¬ 
esting feature that we can see from figure [T51 is that the 
stiffer the EOS the more the quadrupole moment ap¬ 
proaches the Kerr value. This result is well in accor¬ 
dance with previous results that showed that the com¬ 
pactness of the star increases, also the moment of in¬ 
ertia, Love numbers and mass quadrupole approach the 
ones of a black hole, though they will never coincide (see, 
e.g.. Ref. m). Moreover, we confirm in the full rotation 
regime, the previous result obtained in the slow-rotation 
Hartle’s approximation [d^, that the ratio 
is a decreasing function of the NS mass, hence reach¬ 
ing its lowest value at the maximum mass configuration. 
Indeed, as we can see from figure[T31 the largest the max¬ 
imum mass attained by a NS model, the closest the NS 
quadrupole moment approaches the Kerr solution value, 
reaching even values <1.5 for stiff EOS such as the NL3 
model. 


VIII. DISCUSSIONS AND CONCLUSIONS 

We have computed uniformly rotating NSs for se¬ 
lected relativistic mean-field nuclear matter model EOS 
(TMl, GMl, and NL3). Specifically, we have calculated 
their gravitational mass, equatorial and polar radii, ec¬ 
centricity, angular momentum, moment of inertia and 
quadrupole moment. We have established the region 
of stability against mass-shedding and the secular ax- 
isymmetric instability. We have provided plots of all 


these physical quantities e.g. as a function of the mass 
of the configurations. We have also constructed se¬ 
quences of constant rotation frequency and determined 
approximately the rotation rate at which deviations of 
the structure parameters from the spherically symmetric 
(or slowly rotating) values start, obtaining / Ri 200 Hz, 
a value in agreement with previous works (see, e.g., 
Ref. HI). 

From the astrophysical point of view, we have obtained 
a lower bound for the mass of the fastest observed pulsar, 
PSR J1748-2446ad with / = 716 Hz, by constructing its 
constant rotation frequency sequence and constraining it 
to be within the stability region: we obtained Mmin = 
[1.2-1.4] Mq, for the EOS we have used in this work, a 
prediction submitted for observational verification. We 
have obtained also a fitting formula relating the baryonic 
and gravitational mass of non-rotating NSs [see equa¬ 
tion m] independent on the EOS. We have computed a 
formula for the masses of NSs on the secular instability 
line as a function of their angular momentum, see equa¬ 
tion (ED- We studied the Kerr parameter (dimensionless 
angular momentum) of NSs and found that it reaches 
a maximum value (a/M)max ~ 0.7, independent on the 
EOS. This result bring us to the important conclusion 
that the gravitational collapse of a uniformly rotating 
NS, constrained to mass-energy and angular momentum 
conservation, cannot lead to a maximally rotating Kerr 
black hole, which by definition has (a/M)BH.max = 1- We 
have also shown that the quadrupole moment of realistic 
NSs does not reach the Kerr value (for the same val¬ 
ues of mass and angular momentum), but this is closely 
approached from above at the maximum mass value, in 
physical agreement with the no-hair theorem. We have 
also found that the stiffer the EOS the closer the Kerr 
solution is approached. 

It is important to stress that the results shown in 
this work for some specific nuclear EOS likely will re¬ 
main valid in the case of other different models, pro¬ 
viding they are consistent with current observational 
constraints, especially the mass of PSR J0348-I-0432, 
M = 2.O1±O.O4M0 [^. The existence of such a massive 
NS, clearly favors stiff nuclear EOS as the ones obtained 
via RMF theory, which leads to a critical NS mass higher 
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FIG. 13. ratio for the same selected sequences of constant frequency of figure \TI \ and EOS TMl, GMl, and NL3 

(from left to right). We show here only the region of large masses where starts to approach the Kerr value . 


than this constraint. 


To conclude, as already mentioned in sec. IIIIl we would 
like remind the importance of considering a global charge 
neutrality condition for the system, instead of a local 
one, which needs a new and more complete code to treat 
these kind of problems, including the case of fast rotat¬ 
ing strange quark stars with crust, which show similar 
features in the core-crust transition. 
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